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Siiinmarv 

Evolution  of  disturbances  in  a  supersonic  vortex  is  studied  by  using  quasi-cylindrical  linear 
stability  analysis  and  parabolized  stability  equations  (PSE)  formulation.  Appropriate  mean-flow 
profiles  suitable  for  stability  analysis  were  identified  and  modeled  successfully.  Using  the  quasi- 
cylindrical  linear  stability  code,  the  stability  characteristics  of  cylinder  axisymmetric  vortices 
were  mapped  thoroughly  and  discussed  in  detail.  Important  trends  such  as  viscous  and 
compressibility  effects  were  obtained  and  analyzed.  The  results  indicate  that  viscosity  has  very 
little  effect  while  increasing  Mach  number  significantly  stabilizes  the  disturbance.  The 
streamwise  evolution  of  the  constructed  base  flow  was  obtained  by  solving  quasi-cylindrical 
steady-state  equations  using  finite-difference  techniques  and  the  resulting  mean  flow  was 
analyzed  by  using  PSE.  Linear  PSE  analysis  shows  that  the  effect  of  streamwise  mean  flow 
variation  (non-parallelism)  is  small  for  the  case  considered  here.  Nonlinear  evolution  of  helical 
modes  as  well  as  modal  interaction  is  also  studied  by  using  PSE.  It  is  found  that  the  growth  of  the 
disturbances  results  in  the  appearance  of  coherent  large  scale  motion  and  significant  mean  flow 
distortion  in  the  axial  velocity  and  temperature  fields.  In  the  end,  nonlinear  effects  tend  to 
stabilize  the  vortex. 
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1 .  Introduction 

Formation  and  subsequent  downstream  evolution  of  three-dimensional  compressible 
stream  wise  vortices  play  a  significant  role  in  both  external  and  internal  flows.  Effective  control 
or  dissipation  of  vortices  is  an  important  issue  in  such  diverse  areas  as  aircraft  tail  buffeting, 
helicopter  blade  vibration/fatigue,  efficient  mixing  in  combustion  chambers,  and  jet-engine  noise 
reduction.  Unlike  boundary-layer  flows,  however,  the  temporal  and  spatial  evolution  of  com¬ 
pressible  vortices  has  remained  largely  unexplored.  Our  lack  of  understanding  of  the  critical 
features  of  compressible  vortex  evolution,  such  as  stability  characteristics,  laminar-to-turbulent 
transition,  and  compressibility  effects  has  prevented  effective  control  of  vortical  flows  in  practical 
situations. 

Therefore,  the  objective  of  the  present  study  is  to  determine  and  understand  the  stability  char¬ 
acteristics  of  a  physically  realizable  compressible  streamwise  vortex.  This  study  focuses  on 
identifying  the  types  of  instability  modes  present  in  such  a  vortex  with  the  hope  that  these  modes 
may  be  actively  manipulated  in  future  research  on  vortex  control.  In  particular,  the  linear  and 
nonlinear  growth  and  evolution  of  disturbances  in  a  compressible  vortex  is  analyzed  using  quasi- 
cylindrical  linear  stability  analysis  and  parabolized  stability  formulation.  Both  approaches 
complement  each  other  and  together  provide  a  detailed  view  of  the  stability  characteristics  of 
supersonic  vortices. 

2 .  Gk>venmig  Equations 

The  full  governing  equations  for  a  compressible,  variable-property  Newtonian  fluid  flow  in 
cylindrical-polar  (r,0,z)  coordinates  are  given  by  the  Navier-Stokes  equations.  They  are 

mass  conservation: 
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In  equations  (2.1)— (2.5),  Uj.,  Uq,  u^,  p,  T,  and  p  represent  the  three  components  of  velocity, 
pressure,  temperature  and  density,  respectively.  The  coefficients  p,  X,  k  and  Cp  denote  the  fluid 
dynamic  viscosity,  second  coefficient  of  viscosity,  thermal  conductivity,  and  specific  heat  at 


3 


constant  pressure,  respectively,  which  are  taken  to  be  functions  of  temperature  only.  Finally,  the 
above  system  of  equations  is  closed  with  the  addition  of  the  ideal  gas  equation  of  state 

p  =  pRf  .  (2.6) 

3 .  Construction  of  the  Mean-Flow  Profile 

The  selection  of  the  mean-flow  profiles  for  a  supersonic  vortex  is  not  as  straightforward  as  in 
the  incompressible  case  and  therefore  must  be  handled  with  great  care.  One  of  the  main  concerns 
is  to  make  certain  that  the  chosen  base  flow  is  physically  realizable.  Unfortunately,  the  similarity 
solution  of  a  compressible  vortex  obtained  by  Mayer  and  Powell  (1992b)  is  not  a  good  model  for  the 
flow  field  of  interest,  because  their  solutions  do  not  resemble  previous  measurements  of  compress¬ 
ible  vortices.  Therefore,  we  decided  to  construct  the  mean-flow  profiles  from  experimentally 
measured  data.  For  this  purpose,  the  work  of  Smart,  Kalkhoran  and  Bentson  (1994)  provided  some 
of  the  most  extensive  measurements  of  supersonic  streamwise  vortices.  The  construction  of  the 
base  flow  is  presented  and  discussed  in  more  detail  in  the  following  sections. 


3.1  Profile  Based  on  Smart  et  al.  (1994)  Measurements 

The  measurements  of  Smart  et  al.  provide  the  necessary  database  to  construct  the  velocity 
profiles  for  supersonic  wing-tip  vortices.  The  vortex  generation  is  accomplished  via  a  rectangu¬ 
lar  half-wing  with  a  diamond-shaped  cross  section  which  is  mounted  vertically  in  the  test  section. 
Experimental  flow  visualization  revealed  the  supersonic  vortices  to  be  slender.  That  is,  compared 
to  the  radial  gradients,  the  axial  gradients  tend  to  be  small.  Measurements  were  obtained  for 
angles  of  attack  of  10°  (stronger  vortex)  and  5°  (weaker  vortex).  These  vortices  show  significant 
flYifll  velocity  deficit  in  the  vortex  core  and  thus  are  expected  to  be  highly  unstable.  The  10°  case 
was  successfully  modeled  and  is  presented  below. 

The  procedure  for  constructing  the  profiles  is  as  follows.  In  the  cylindrical  coordinates 
let  n(r)  and  7}(r)  represent  the  non-dimensional  pressure  and  temperature,  respectively.  In 
what  follows,  the  pressure,  temperature,  and  density  are  normalized  by  the  far-field  values  of 
d5mamic  pressure  ,  static  temperature  and  density  p„,  respectively.  The  radial  coordi¬ 

nate  is  normalized  by  the  value  of  the  radial  coordinate  Tq  where  the  tangential  velocity  is  a 
maximum,  and  the  velocities  are  nondimensionalized  by  the  far-field  value  of  the  axial  velocity 
Woo. 

Noting  that  for  a  slender  vortex  the  radial  velocity  is  negligible,  good  fits  to  the  following 
expressions  for  the  measured  tangential  and  axial  velocities  were  foimd  to  fit  the  data  well. 


Ue=V{r)  =  ^\l-e-'-  J, 
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and 


u,  =  W{r)  =  \-5e-'^-^^\ 


(3.2) 


where  q  is  the  usual  swirl  parameter  similar  to  that  defined  for  the  Batchelor  vortex  (see  Lessen  et 
al.  1974)  and  5  is  the  axial  velocity  deficit.  Unfortunately,  no  measurements  of  temperature 
across  the  vortex  were  reported  by  Smart  et  al.  (1994).  As  a  first  step,  the  total  temperature  is 
assumed  to  be  constant  throughout  the  vortex.  The  total  temperature  1]^  then  is  given  by 

T](  =  constant  =  p  (3.3) 


where  at  a  sufficiently  large  radius,  P  is  determined  from 


P  =  -l  +  i{Y-l)Ml 


with 


(3.4) 


Assuming  the  flow  is  locally  adiabatic,  the  static  temperature  distribution,  rf  ,across  the  vortex  is 
given  by 


rj{r]  =  p-^^Ml(v^  +  W^). 


(3.5) 


For  a  quasi-cylindrical  vortex,  the  radial  momentum  equation  for  the  mean  flow  reduces  to 


(3.6) 
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Integrating  both  sides,  we  obtain 


n(r)  = 
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A  numerical  integration  thus  provides  the  pressure  after  substituting  the  experimentally 
determined  Mq  distribution  in  equation  (3.8).  Finally,  the  density  variation  is  obtained  from  the 

ideal  gas  law 

.  (3.9) 

ri 

3.2  Quasi-Cvlindrical  Computation  of  Base  Flow 

The  fitted  vortex  (3.1)-(3.9)  is  used  as  an  initial  profile  to  obtain  the  streamwise  evolution  and 
quashcylindrical  nature  of  the  vortex  base-flow.  The  computed  mean  flows  at  different 
streamwise  locations  were  subsequently  used  to  determine  the  effect  of  core  expansion  on  the 
vortex  stability  characteristics.  In  addition,  the  computed  mean  flows  were  used  in  conjunction 
with  the  Parabolized  Stability  Equation  (PSE)  formulation  to  determine  linear  and  nonlinear 
evolution  of  disturbances  in  supersonic  vortices. 

A  numerical  method  for  modeling  isolated  compressible  vortices  is  presented.  The  formula¬ 
tion  is  based  on  the  assumptions  that  the  vortex  is  steady,  quasi-axisymmetric,  and  slender.  The 
governing  equations  and  the  solution  procedure  are  discussed  in  the  following  sections. 

3.2.1  Governing  equations 

The  governing  equations  are  derived  starting  from  the  full  unsteady,  compressible  Navier- 
Stokes  equations  in  cylindrical  coordinates  (given  by  Eqs.  (2.1)-(2.5)  in  Section  2.0)  where  z  is  the 
axial  direction,  r  is  the  radial  direction,  0  is  the  azimuthal  direction,  and  u,  v,  w  are  the  velocity 
components  in  the  r,  0,  and  z  directions,  respectively.  Assuming  the  vortex  is  steady  (3/3t  =  0), 
quasi-axisymmetric  (0/00  =  0),  and  slender  (fcore  «  L  where  Fcore  is  the  vortex  core  radius  and  L 
is  the  breakdown  length),  the  Navier-Stokes  equations  reduce  to 

Continuity: 

Axial  momentum: 

Azimuthal  momentum: 

Radial  momentum: 
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Energy; 


Equation  of  state: 


pCp 


8T  8T 
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p  =  pRT 


(3.15) 


Here,  p  is  the  density,  p  is  the  static  pressure,  T  is  the  static  temperature,  R  is  the  gas  constant,  and 
|1,  k,  and  Cp  are  the  coefficients  of  dynamic  viscosity,  thermal  conductivity,  and  specific  heat  at 
constant  pressure,  respectively.  This  set  of  equations  is  similar  to  the  boundary-layer  equations 
except  that  the  radial  momentum  equation  is  retained  and  the  pressure  is  allowed  to  vary  through 
the  vortex  core. 

Equations  (3.10)  through  (3.15)  are  nondimensionalized  based  on  the  following  procedure: 
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2l. 

fo 
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T 

ai/Cp 


Re  = 


where  Tq  is  the  initial  vortex  core  radius,  is  the  freestream  speed  of  sound,  p«,  is  the  fireestream 
density,  is  the  freestream  coefficient  of  dynamic  viscosity,  is  the  freestream  Mach  number, 
Pr  is  the  Prandtl  number,  and  Re  is  the  Re3molds  number.  In  terms  of  these  nondimensional 
variables,  the  governing  equations  become  (the  overbars  are  omitted) 


Continuity: 
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Azimuthal  momentum: 
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Radial  momentum: 
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Energy: 
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Equation  of  state: 


dT  ai 
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In  the  present  formulation,  \l  is  determined  from  Sutherland's  viscosity  law  which  takes  the 
following  form  in  terms  of  the  nondimensional  variables: 

H  =  [(Y  -  DTf 


Once  |J.  is  known,  the  thermal  conductivity  k  is  determined  from  the  Prandtl  nvunber.  The  bound¬ 
ary  conditions  at  r  =  0  and  r  — >  <»  are 


r  =  0:  u  = 


(3.23) 


r  — >  oo: 


dw  _  d(rv)  _  dT 
dr  dr  dr 


(3.24) 


p  =  Pe(z)  (3.25) 

With  M„„,  Re,  T„,  Pr,  and  Pe(z)  specified,  Eqs.  (3.16)-(3.21)  can  be  marched  along  the  axial 
direction  starting  from  a  set  of  initial  profiles  for  w,  v,  and  T  at  z  =  Zq.  The  initial  profiles  for  p 
and  p  are  obtained  from  Eqs.  (3.19)  and  (3.21)  while  u  is  set  to  zero  at  z  =  Zq. 

3.2.2  Solution  procedure 

Equations  (3.16)-(3.21)  are  solved  in  a  decoupled  manner  at  each  marching  step.  First,  the 
axial  momentum,  azimuthal  momentum,  and  energy  equations  are  marched  to  obtain  w,  v,  and 
T.  Next,  p  and  p  are  computed  by  integrating  the  radial  momentum  equation  together  with  the 
equation  of  state.  Finally,  u  is  obtained  by  integrating  the  continuity  equation.  The  entire 
solution  procedure  is  described  in  detail  below. 
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step  1: 

In  the  first  step  w,  v,  and  T  are  computed  from  Eqs.  (3.17),  (3.18),  and  (3.20).  These  three 
equations  are  first  written  in  vector  form  as  follows 
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Equation  (3.26)  is  discretized  using  the  implicit  Crank-Nicolson  scheme  which  is  second-order 
accurate  in  both  spatial  directions.  On  a  computational  grid  with  uniform  grid  spacing,  the 
resulting  expression  is 
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Here,  superscript  i  and  subscript  j  denote  the  discrete  axial  and  radial  locations,  respectively,  and 
superscript  n  represents  the  nonlinear  iteration  index.  The  solution  of  this  equation  requires  the 
inversion  of  a  3  x  3  block  tridiagonal  matrix.  The  boundary  conditions  applied  at  r  =  0  and  r  =  Tmax 
are 

r  =  0:  —  =  0  — >  Wi  —  W2 

dr 

V  =  0  ^  Vi  =  0 
^=0^Ti=T2 
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Step  2: 

After  w,  V  and  T  are  determined,  p  is  computed  by  integrating  Eq.  (3.19)  together  with  Eq. 
(3.21)  from  r  =  Tmax  to  r  =  0.  At  r  =  Fmax.  P  is  obtained  from  the  equation  of  state  as  follows; 

D 

Pjmax-  .  y 

/  “  i  A  jxnax 

Here  either  Pe  or  dpjdz  needs  to  be  specified  at  r  =  fmax-  If  dpjdz  is  specified,  Pe  is  obtained  as 


pie+1  =  pi  +  Az 


ape 

dz 


Once  p  is  known,  p  is  computed  from  the  equation  of  state. 


Step  3: 

In  the  final  step,  u  is  computed  by  integrating  Eq.  (3.16)  from  r  =  0  where  u  =  0  to  r  =  fmax  using 
the  values  of  w  and  p  obtained  in  steps  1  and  2,  respectively. 

After  p,  p,  and  u  are  computed  as  described  in  steps  2  and  3,  the  left  hand  side  matrix  A  and  the 
right  hand  side  vector  B  in  Eq.  (3.27)  are  then  updated  and  the  entire  process  (i.e.,  steps  1  through  3) 
is  repeated  until  the  solution  converges.  Note  that  in  this  solution  procedure,  the  streamwise 
pressure  gradient  term  in  Eq.  (3.26)  is  specified  and  treated  as  a  source  term  at  each  iteration. 
Therefore,  the  stability  problem  associated  with  space  marching  in  subsonic  flow  regions  is 
avoided. 

Numerical  computations  of  the  mean  flowfield  allows  the  variation  of  parameters  such  as 
vortex  strength  q,  Mach  number  M^e,  euid  Reynolds  number  Re  independently.  In  this  manner, 
important  effects  (e.g.,  compressibility,  viscous,  etc.)  can  be  isolated.  Smart  et  al.’s  experiments 
were  conducted  at  a  nominal  Mach  number  of  2.49.  Although  variation  of  different  flow 
parameters  may  cause  the  computed  flowfield  to  become  unrealizable  in  their  tunnel,  the 
realizability  of  the  flowfield  in  some  other  tunnels  is  insured  as  long  as  numerically  converged 
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solutions  are  obtained.  The  computed  profiles  for  q  =  0.2,  Re  -  30000,  and  Moo  -  2.49  at  several  axial 
locations  are  presented  in  Figures  la-lg.  The  vortex  initial  axial  velocity  deficit  is  35%  for  these 
cases.  Note  viscous  diffusion  only  significantly  affects  axial  velocity,  static  temperature,  and 
density  while  tangential  velocity  and  static  pressure  remain  unaffected. 

4.  Stability  Equations 

The  set  of  equations  (2.1)  -  (2.6)  provide  the  starting  point  for  the  derivation  of  any  form  of 
stability  equations.  In  the  present  case,  the  base  flowfield  obtained  by  solving  the  axisymmetric 
slender  vortex  equations  (section  3.2)  is  perturbed  by  small  fluctuations.  The  flow  variables  then 
are  assumed  to  consist  of  a  mean  part  and  a  fluctuating  part,  i.e., 

u  =  U^+eUr,  v  =  Ug+eug,  w  =  U^+eu'^ 

p  =  P  +  sp\  p^p  +  ep',  f  =  T  +  ^'  (4.1) 

and 

p  =  p(f),  A  =  A(f),  r  =  (4.2) 

where  e  denotes  an  infinitesimal  quantity.  Substituting  equation  (4.1)-(4.2)  into  the  (2.1)-(2.6) 
and  collecting  terms  in  powers  of  e,  one  obtains  the  equations  governing  linear  and  nonlinear 
evolution  of  small  disturbances  in  compressible  streamwise  vortices.  These,  equations  which  are 
quite  complicated  and  long  (but  straight  forward  to  derive),  provide  the  basis  for  developing  the 
linear  stability  and  parabolized  stability  formulations. 

4.1  Qiiasi-Cvlindrical  Linear  Stability  Formulation 

Linear  stability  analysis  is  applicable  in  the  early  stages  of  the  disturbance  evolution  where, 
due  to  the  smallness  of  the  fluctuation  amplitude,  no  significant  modal  interaction  or  mean  flow 
distortion  exists.  The  governing  equations  are  obtained  by  keeping  terms  of  0(e)  and  neglecting 
all  higher  order  terms.  In  this  formulation,  the  base  flow  is  very  slowly  developing  in  the 
streamwise  direction  and  therefore  at  each  downstream  z  station  the  local  mean  flow  quantities 
are  taken  to  be  functions  of  the  radial  coordinate  only.  Next,  the  perturbations  are  assumed  to  be 
composed  of  three-dimensional  helical  Fourier  modes.  The  flow  variables  are  then  written  as 
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^i{az-\-n6-(at) 


(4.3) 


u 

mr)' 

Fir)' 

V 

V{r) 

G(r) 

w 

W(r) 

H(r) 

p 

n{r) 

Pir) 

f 

nir) 

Tir) 

p. 

pir) 

Air) 

where,  F,  G,  H,  P,  T  and  A  are  the  disturbance  complex  amplitudes  which  are  only  functions  of  the 
radial  coordinate.  For  the  temporal  solution,  the  variables  a  and  n  are  the  real-valued  wavenum¬ 
bers  in  the  axial  and  azimuthal  directions,  respectively.  The  azimuthal  wave  number  is  an 
integer  where  n  9^  0  indicates  helical  waves  and  n  =  0  represent  axisymmetric  modes.  The 
angular  frequency  (O  is  complex,  and  the  sign  of  its  imaginary  part,  ©j,  determines  whether  the 
disturbances  will  grow  (imstable  if  (Oi  >  0)  or  decay  (stable  if  Wi  <  0).  For  the  spatial  formulation, 
CO  and  n  are  taken  to  be  real  and  complex  a  is  sought.  If  the  imaginary  part  of  the  wave  number  a 
is  negative  ai  <  0,  the  disturbance  grows  in  space  (unstable),  and  if  aj  >  0  it  decays  (stable). 
Substituting  (4.3)  into  the  Oie)  terms  in  the  equations  and  neglecting  terms  containing 
streamwise  derivative  of  the  mean  flow,  the  linearized  stability  equations  in  terms  of  the 
perturbation  eigenfunctions  are  obtained.  These  quasi-cylindrical  equations  in  nondimensional 
form  are: 

mass  conservation: 


dr 


1  dn 
n  dr 


1  dr]  ^ 

T}  dr  r 


lint  r*  Irr  f  G 1  rfP 
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(4.4) 


r-momentum: 
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\d^Xdr]fdU  U\^\dl(d}U  U  ^ldU\^indii(\dV  V\ 
^  fi  dif'  rfr  V  rfr  r  J  fi  dt]  ^  dr^  r  r  dr  ^  ^  dr) \r  dr  r  ) 

+^fl7^_ZlV— +  -^f-— -— ll7’  =  0 

fir)  dr  r  j  n  dr)  dr  n  dr)\r  dr  r^JJ 


^-momentum: 


RepdV  RepV  ^  in  I  dp  drj  ^  ^  d^G  ^  1  dja  dr]  RepU 
p  dr  pr  r  p  dr]  dr  ^  dr^  _p  dr]  dr  p 


lldG  icoRep  inRepV  iaRepW  Repf7  ^ 2  \  ^  ^  dp  dr]  1  ^ 

Ji  fi  Itr  r^[  r  ^  dr)  dr  r^ \ 
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^  \  dpU'^V  IdV^^  2  dpfdV  ^"[L-o 
pd^ydr'^  r^  r  dr  j  rpdpl^rfr  rJJ 


(4.6) 
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^-m  omentum: 
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energy-equation: 
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T  =  0 


(4.8) 


Ky-mlVr  duUdU^  U^\ 

K  di\  y  dr )  j 


^  (y-DMlPrdu  \(dV  ^  (y-l)MlPr  dX  \(dU  ^  t/^] 

K  dr]  \  dr  r )  ^  k  dv]  \  dr  r  j  ^ 


In  the  linearized  governing  equations,  the  parameters  Re,  Pr  and  represent  the  flow  Reynolds, 
Prandtl,  and  Mach  numbers,  respectively.  Moreover,  fluctuations  in  the  flow  properties  are 
related  to  the  temperature  perturbations  via  the  assumed  functional  form  of  viscosity  and  thermal 
conductivity. 


4.1.1  Boundary  Conditions  and  the  Numerical  Method 

The  compressible  governing  equations  (4.4)-(4.8)  require  eight  boundary  conditions.  The 
appropriate  boundary  conditions  on  the  centerline  of  the  vortex  depend  on  the  wave  number  n,  and 
must  be  chosen  such  that  the  solutions  remain  single-valued  and  smooth.  The  conditions  at  r  =  0 
(Eliorrami  1991a)  are 

if  n  =  0,  F(0)  =  G(0)  =  0 


if  71  =  ±1,  F(0)  +  7G(0)  =  0 


ff(o)=r(o)=o 


(4.9) 


if  |71|>1,  i?’(0)  =  G(0)  =  0 

2?(0)  =  T(0)  =  0. 

The  other  set  of  the  boundary  conditions  are  applied  at  a  large  radial  distance  from  the  vortex 
centerline.  Normally,  at  the  far-field  boundary,  the  disturbances  are  taiken  to  approach  zero;  that 
is, 

asr-^oo  F,G,H,T,^Q  .  (4.10) 


The  proper  far-field  behavior  of  the  disturbances  is  shown  more  formally  below.  We  restrict  our 
attention  to  an  inviscid  analysis,  even  though  in  the  numerical  computations  the  effect  of  viscosity 
is  fully  taken  into  account.  The  reason  for  this  restriction  is  two-fold.  First,  from  the  previous 
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incompressible  results  for  three-dimensional  vortices,  the  inviscid  disturbances  are  known  to  be 
more  unstable  than  the  purely  viscous  modes  of  instability  (Lessen  et  al.  1974,  Mayer  and  Powell 
1992a,  Khorrami  1991b).  Second,  the  Re5uiolds  numbers  of  the  experiments  are  quite  high  and 
therefore  the  inviscid  disturbances  are  expected  to  dominate. 

The  limiting  form  of  the  mean-flow  field  developed  in  section  3.1  as  r  ->  <»  is 


W,  7],  p->l. 


77 


r 


dV  q 
dr  ^ 


(4.11) 


dW  ^  ^  dn  ^  Q 

dr  ’  dr'  dr’  dr 

Substituting  these  quantities  into  the  governing  stability  equations  for  an  inviscid,  non-heat¬ 
conducting  gas,  the  as5unptotic  form  of  the  governing  equations  in  the  far-field  are  obtained.  After 
much  simplification,  we  get 


d^H 

1 

fl  q^Ml' 

dr^ 

z' 

^  1 
1 

]dr  [ 

2q^Ml 


77  =  0  . 


J 


(4.12) 


At  a  large  radial  distance  from  the  vortex  centerline,  the  solution  to  (4.12)  approaches 


(4.13) 


where  is  the  modified  Bessel  fiinction  of  the  second  kind  of  order  n,  and  a  is  given  by 


<7  = 


a?'  - Ml,[a - (of 


1/2 


(4.14) 


The  real  part  of  o  is  taken  to  be  positive  for  the  solution  to  remain  finite.  From  equations  (4.12)- 
(4.14),  notice  that  in  the  limit  as  >  0 ,  the  incompressible  solutions  (Lessen  et  al.  1974,  Ash  and 
Khorrami  1995)  are  recovered.  Proceeding  a  step  further,  it  can  be  shown  that  for  large  arguments, 
the  leading  term  expansion  in  (4.13)  is 

•  (4.15) 

or 

The  above  behavior  clearly  demonstrates  that  the  disturbances  have  a  faster  than  exponential 
decay  rate  as  r  ^  oo.  Thus,  we  may  truncate  the  computational  domain  at  a  large  radius 
enforce  the  boundary  conditions  given  by  (4.10). 
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The  equations  (4.4M4.8)  and  the  boundary  conditions  (4.9M4.10)  constitute  a  well-posed 
eigenvalue  problem.  Notice  that  there  are  no  pressure  boundary  conditions  given  because  the 
governing  equations  are  discretized  using  a  Chebyshev  spectral  collocation  technique  imple¬ 
mented  on  a  staggered  grid.  In  this  method,  the  pressure  terms  are  evaluated  only  at  the  interior 
nodal  points,  thus  eliminating  the  need  for  any  pressure  conditions.  The  method  is  an  extension 
of  the  formulations  of  Macaraeg,  Streett  and  Hussaini  (1988)  and  Khorrami  (1991c)  to  the 
compressible  rotating  flow  regime.  The  technique  is  very  efficient  and  has  been  successfully 
applied  to  the  case  of  compressible  swirling  jets  by  Khorrami  (1991a,  1995).  Next,  for  the  temporal 
formulation,  the  discretized  equations  are  arranged  in  a  generalized  eigenvalue  format 

£jX  =  ft)E2X,  (4.16) 

where  Ej  and  are  the  square  coefficient  matrices  of  and  Np  is  the  number  of  Chebyshev 

polynomials  used  in  the  expansion  of  the  eigenfunctions.  The  eigenvector  x  is 

x  =  [F  G  H  P  rf  ,  (4.17) 

where  the  superscript  T  denotes  the  transpose.  The  solution  of  (4.16)  for  the  eigenvalue  co  is 
obtained  utilizing  any  readily  available  QZ  or  LZ  routines. 

For  the  present  study,  the  physical  coordinate  [0,rjj^ax]  is  mapped  into  the  computational 
domain  [-1,1]  using  the  algebraic  transformation 

r  =  ai^  ,  (4.18) 

b-^ 

where  a  and  b  are  constants  with 

^  =  1  +  ^^.  (4.19) 

^ax 

Throughout  our  studies,  we  found 

4  <  a  <  10  ,  (4.20) 

and 

20<rniax^60,  (4.21) 
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to  be  sufficient  to  ensure  that  the  computed  results  are  domain  independent  and  well  resolved. 
Furthermore,  the  number  of  collocation  points  needed  to  obtain  six  or  seven  significant  digits  of 
accuracy  varied  between  70  <  Np  <130  , 

4.2  Parabolized  Stability  Formulation 

Parabolized  stability  formulation  was  pioneered  by  Herbert  (1991),  Bertolotti  (  1991),  Bertolotti 
et  al.  (1992),  Chang  et  al.  (1991),  £uid  Chang  and  Malik  (1994).  The  approach  is  applicable  to  flow 
fields  that  meet  the  following  two  important  criteria.  First,  the  base  flow  variations  in  the 
streamwise  direction  are  relatively  small  and  second,  the  disturbances  are  convective  in  nature. 
For  supersonic  streamwise  vortices,  which  are  of  interest  here,  both  of  the  above  conditions  are 
met. 


4.2.1  Governing  Equations 

The  governing  equations  for  PSE  are  obtained  by  keeping  all  terms  in  the  equations  (A.l)  •“ 
(A,5)  (appendix  A)  up  to  the  0(£^).  The  full  nonlinear  equations  are  grouped  and  written  in 
general  form  as 


.90  ^90  C 

12  — +  A— h-B— +  — +  Z><P 
9t  9z  9r  r 


9^0 

9z^ 


9^0 


Q.b 

r  9z9G 


Q^g  d^0  Qgg  d^<p' 

r  9r90  96^ 


(4.22) 


where  0  is  the  perturbation  vector  given  as 
^  =  {p'.Wz,Wr.“0>^'}^. 

£uid  Re  is  the  Reynolds  number  as  defined  in  section  4.1.  Matrices  C,  D, 

Qj.q,  and  Qqq  are  Jacobians  of  the  corresponding  total  flux  vectors  which  are  composed  of  a  linear 
part  with  only  mean  flow  quantities  (denoted  by  superscripts  £  below)  and  a  nonlinear  part  which 
contains  perturbation  quantities  (denoted  by  superscripts  n£  below). 

As  in  the  previous  section,  the  perturbations  are  assumed  to  be  periodic  in  time  and  azimuthal 
direction  resulting  in  an  expansion  of  the  form 


Hr,9,z)=  ^'^Xmn{r,z)e' 


i(n9~mot) 


(4.23) 
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For  most  flows  of  interest,  it  is  sufficient  to  truncate  the  Fourier  decomposition  (4.23)  after  a  finite 
number  of  modes 

M  N 

4.(r,9,*)=  ^  (4.24) 

m=~M  n=-N 


where  M  and  N  are  half  the  number  of  modes  kept  in  the  truncated  series. 

We  are  interested  in  the  evolution  of  wave-like  disturbances  that  are  growing  or  decaying 
while  convected  downstream.  Thus  it  can  be  assumed  that  is  composed  of  a  fast  oscillatory 

wave  part  plus  a  slowly  var3dng  shape  function,  i.e., 


Xmn{r,z)  =  Wmn  {r,z)E^n{z), 


(4.25) 


E^mn{z)- 


(4.26) 


Substituting  (4.24)-(4.26)  into  the  governing  equations,  we  obtain  the  equations  for  the  shape 
function  as 


^mnW mn  -^mn 


dyf, 


mn  D  _ 

dz  dr 


Re 


Qi 


dzdr 


dr^ 


=  F^JE^,.  (4.27) 


The  left  hand  side  of  equation  (4.27)  contains  only  linear  coefficient  matrices  and  all  nonlinear 
terms  are  included  in  the  forcing  function  which  is  the  Fourier  component  of  the  total  forcing, 
defined  as 


where  all  matrices  with  superscript  nt  denote  corresponding  coefficients  containing  dependent 
variables  0  (nonlinear  terms).  The  forcing  function  is  defined  in  physical  space  and  can  be 
evaluated  by  a  similar  Fourier  series  expansion 
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(4.29) 


M  N 


F'^^{r,e,z)= 


9-m(Ot) 


m=-M  n-~N 


Computationally,  the  Fourier  decomposition  of  equation  (4.29)  can  be  accomplished  by  using  the 
Fast  Fourier  Transform  (FFT)  of  F"'  which  is  evaluated  numerically  in  the  physical  space.  The 
corresponding  linear  formulation  is  obtained  by  setting  F^n  to  zero  in  (4.27). 

As  it  stands,  equation  (4.27)  is  elliptic  in  nature  in  particular  due  to  the  presence  of  the  term 
— The  computational  cost  of  solving  such  an  equation  is  prohibitively  expensive.  PSE 
formulation  removes  this  source  of  ellipticity  based  on  physical  and  dimensional  arguments.  A 
brief  overview  of  the  procedure  is  given  below. 

In  the  PSE  formulation,  the  perturbation  quantities  for  a  particular  mode  are  expanded  as, 


.  \  i\  (xdz-^inQ-icot 

u^=G(r,z)e^ 


(4.30) 


with  similar  expansions  for  other  flow  variables.  Here  c.c.  represents  complex  conjugate  to 
ensure  the  physical  quantities  are  real  and  is  dropped  from  the  subsequent  discussion.  Note  that 
unlike  the  quasi-cylindrical  expansion  (4.3),  here  the  disturbance  amplitude  is  dependent  on  both 
radial  and  axial  coordinates  with  the  underlying  assumption  that  streamwise  variation  is  slow. 
Constructing  the  first  and  second  derivatives  of  G  in  z,  we  have 


(4.31) 


-f  =  ya^(z)G(r,z)  + 


.daiz)  /  v(9G  d^G]  i\cdz-^ine- 


(4.32) 


Since  the  mean  flow  variations  are  0(Re  M,  for  moderate  and  high  Re3molds  number  flows  the 

d^G 

wave  amplitude  is  weakly  dependent  on  the  z  coordinate.  For  these  flows  the  term  — ^  cam  be 

oz 

omitted  amd  the  second  derivative  is  approximated  as 


^  =  \-a^z)G{r,z)  . 


(4.33) 


Substituting  expansion  (4.30)  into  the  equation  (4.27)  and  utilizing  approximation  (4.33),  all  the 
second  derivative  terms  will  be  removed  and  a  set  of  partial  differential  equations  which  contain 
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only  first  derivatives  in  z  is  obtained.  With  the  above  simplification  and  after  neglecting  all 
terms  of  O^Re”^  j,  the  following  reduced  equation  is  obtained 


d¥mn 

dz 


^¥mn 

dr 


Qrr  ^  ¥ mn  _  p 

Re  dr^ 


/E^ 


(4.34) 


The  emerging  equation  (4.34)  which  describes  the  evolution  of  the  shape  fimction  is  almost 
parabolized.  For  complete  parabolizing  (4.34),  one  needs  to  address  the  convective  side  of  the 
governing  equation.  Here,  the  elliptic  effects  associated  with  the  wave  part  is  absorbed  in  matrices 
Dmj^,  and  The  remaining  ellipticity  in  the  governing  equation  caused  by  the  terms  — 

is  a  non-issue  for  the  supersonic  vortices  of  interest  due  to  the  absence  of  a  subsonic  region  in  the 
flow  field. 

4.2.2  Boundary  Conditions 

The  boxindary  conditions  on  the  centerline  r  =  0,  are  identical  to  the  conditions  (4.9)  derived 
for  the  quasi-cylindrical  stability  formulation.  In  the  farfield  as  the  Dirichlet  boundary 

conditions  (4.10)  can  be  employed.  For  supersonic-mode  disturbances,  however,  the  farfield 
Dirichlet  conditions  may  not  be  appropriate  since  the  disturbance  does  not  decay  rapidly  outside  of 
the  vortex  core.  Therefore,  for  both  linear  and  nonlinear  PSE  calculations,  the  non-reflecting 
boxmdary  conditions  of  Thompson  (1987)  are  applied  in  the  farfield. 

4.2.3  Numerical  Approach 

Numerical  solution  of  the  parabolized  stability  equation  (4.34)  requires  discretization  in  both 
r  and  z  directions.  The  discretization  procedure  is  identical  to  the  approach  employed  by  Malik 
and  Chang  (1997)  in  their  PSE  study  of  a  supersonic  axisymmetric  jet  flow.  That  paper  describes 
the  intricacies  of  the  numerical  method  in  detail  and  therefore  is  not  repeated  here.  For  the  present 
study  the  farfield  boundary  conditions  were  applied  at  r^^  =  25  which  proved  to  be  more  than 
sufficient. 

5 .  Results  and  Discussions 

The  governing  equations  are  written  in  a  general  form  valid  for  a  perfect  gas  with  variable 
properties.  The  temperature  dependency  of  the  fluid  properties  may  be  assumed  to  follow 
Sutherland’s  law  or  a  simple  power  law  variation.  In  the  present  work,  for  the  quasi-cylindrical 
linear  stability  calculations,  we  obtained  results  with  both  constant  and  variable  flow  properties. 
We  found  the  difference  between  the  two  sets  of  computations  (in  terms  of  the  computed  growth  rate 
and  frequency)  to  be  less  than  3%.  Therefore,  for  quasi-cylindrical  linear  stability  analysis  only 
results  obtained  with  constant  flow  properties  are  presented.  For  both  linear  and  nonlinear  PSE 
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computations  Sutherland’s  law  was  assumed  with  a  freestream  static  temperature  of  130°K. 
Overall,  for  the  Mach  number  range  of  interest  these  variations  do  not  play  a  significant  role.  In 
addition,  the  second  coefficient  of  viscosity  is  taken  to  be  A  =  -2/3jU. 

The  parameter  space  affecting  stability  is  multidimensional.  We  expect  the  stability  charac¬ 
teristics  of  supersonic  vortices  to  be  dependent  on  g,  5,  Af,  Re,  Pr,  a,  (0,n,  etc.  For  the  problem  to  be 
tractable,  some  restriction  on  the  values  of  these  parameters  is  warranted.  The  Prandtl  number  Pr 
was  set  to  0.76  to  match  the  Smart  et.  al  experiment.  The  Reynolds  number  Re  was  chosen  to  be  as 
in  the  Smart  et.  al  experiment.  Owing  to  the  inviscid  nature  of  the  disturbance  field,  for  moderate 
and  high  Re,  the  Reynolds  number  does  not  play  a  significant  role.  The  vortex  initial  axial 
velocity  deficit  5  is  35%  which  matches  the  experiment.  As  in  the  case  of  incompressible  vortices 
(see  Lessen  et.  al  1974),  the  ratio  q!  8  determines  stability  or  instability  and  therefore  not  much 
generality  is  lost  by  fixing  5. 

Numerical  computations  of  the  mean  flowfield  allow  one  to  vary  parameters  such  as  vortex 
strength  q,  Mach  number  Moo,  and  Reynolds  number  Re  independently.  In  this  manner, 
important  effects  (e.g.,  compressibility,  viscous,  etc.)  can  be  isolated.  Smart  et  al.’s  experiments 
were  conducted  at  a  nominal  Mach  number  of  2.49.  Although  variation  of  different  flow 
parameters  may  result  in  a  computed  flowfield  that  is  unattainable  in  their  tunnel,  the 
realizability  of  the  flowfield  in  some  other  tunnels  is  insured  as  long  as  numerically  converged 
solutions  are  obtained.  The  computed  profiles  for  q  =  0.2,  Re  =  30000,  and  Moo  =  2.49  at  several  axial 
locations  are  presented  in  Figures  la-lg.  The  vortex  initial  axial  velocity  deficit  is  35%  for  these 
cases.  Note  viscous  diffusion  only  significantly  affects  axial  velocity,  static  temperature,  and 
density  while  tangential  velocity  and  static  pressure  remain  unaffected.  Table  I  presents  values 
and  ranges  of  different  parameters  used  in  the  present  study. 

5.1  Quasi- Cylindrical  Linear  Stability  Analysis 

The  local  linear  stability  analysis  is  very  adept  at  determining  the  relevant  parameter  space, 
general  trends,  and  overall  stability  characteristics.  These  computations  also  provide  the  guid¬ 
ance  to  conduct  the  nonlinear  PSE  calculations.  We  now  proceed  with  presenting  the  computed 
results.  The  reference  point  to  be  established  for  this  vortex  is  the  case  of  g  =  0.  Here  the  flow  is 
simply  that  of  an  axisymmetric  wake  which  is  only  unstable  with  respect  to  instabilities  having 
an  azimuthal  wave  number  of  zi  =  ±1 . 

The  variation  of  growth  rate  versus  frequency  of  the  n  -  1  mode  for  the  wake  is  displayed  in 
figure  2.  The  instability  has  a  modest  growth  rate  which  does  not  change  significantly  as  the  wake 
evolves  downstream.  The  growth  rate  curve  has  a  peak  value  of  0.027  at  the  initial  location  which 
subsequently  drops  to  0.022  at  a  station  320  core  radii  downstream.  The  nondimensional 
frequency  co^  corresponding  to  is  0.8.  In  contrast,  the  vortex  with  a  small  amount  of  swirl  (g 
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=  0.1)  is  highly  destabilized  by  the  asymmetric  /i  =  1  disturbance  as  shown  in  figure  3. 
Specifically,  the  growth  rate  curves  depict  a  tripling  effect  when  compared  with  the  wake 
instability.  In  addition,  there  is  nearly  a  20%  reduction  in  the  instability s  growth  rate  at  z  =  320 
compared  to  the  initial  station  z  =  0.  The  effect  of  increased  swirl  on  the  n  =  1  mode  is  displayed  in 
figures  4—6.  The  instabilitys  amplification  factor  peaks  at  g  =  0.2  and  there  is  a  significant 
reduction  in  <Wj  as  the  disturbance  travels  downstream.  Increasing  q  beyond  0.2  has  a  stabilizing 
influence  on  the  n  =  1  mode.  The  growth  rate  for  q  =  0.4  is  comparable  to  the  values  obtained  for  the 
wake  disturbances  (figure  2).  At  this  q  the  n  =  1  mode  is  unstable  over  a  narrow  band  of 
frequencies.  Moreover,  the  perturbation  is  stabilized  rapidly  as  it  moves  in  the  axial  direction  and 
at  a  z  location  slightly  greater  than  70  becomes  a  stable  mode. 

Unlike  the  axisymmetric  wake  flow,  the  vortex  is  unstable  with  respect  to  higher  azimuthal 
wave  numbers.  The  amplification  behavior  of  the  n  —  2  mode  for  values  ofq  =  0.1,  0.2,  0.3,  and  0.4 
are  shown  in  figures  7-10.  Atq  =  0.1,  the  instability  is  growing  at  a  lower  rate  than  the  n  =  1  distur¬ 
bances  (figure  7).  At  q  =  0.2,  the  growth  rates  for  this  mode  become  comparable  to  the  n  =  1  distur¬ 
bance  and  are  slightly  lower  only  in  the  first  few  axial  locations.  On  the  other  hand,  the  growth 
rate  of  the  n  =  2  mode  virtually  remains  unchanged  for  large  values  of  z.  This  behavior  becomes 
more  pronounced  at  higher  swirl  values.  Thus,  while  the  n  =  1  disturbances  may  dominate  the 
initial  stages  of  the  Unear  growth,  it  is  the  n  =  2  modes  that  tend  to  dominate  at  farther  locations 
downstream.  Similar  to  the  n  =  1  perturbations,  the  growth  rate  for  the  n  =  2  mode  peaks  atq  =  0.2. 
Further  increases  in  the  swirl  parameter  provide  a  stabilizing  effect  up  to  a  q  =  0.5  where  the 
vortex  becomes  totaUy  stabiUzed  with  respect  to  this  perturbation. 

The  vortex  is  also  unstable  with  respect  to  larger  azimuthal  wavenumbers.  We  have  mapped 
thoroughly  the  stability  characteristics  of  the  vortex  for  n  =  3,  and  4.  Figure  11  displays  the 
variation  of  coi  vs  the  nondimensional  frequency  co^  for  the  n  =  3  mode  at  q  =  0.2.  Unlike  the 
previous  two  modes,  the  maximum  growth  rate  does  not  occur  at  the  initial  station.  Rather,  the 
instability  attains  its  peak  amplification  rate  at  a  downstream  location  of  z  =  120.  Beyond  this 
axial  station  the  amplification  rate  diminishes  at  a  very  slow  pace  as  the  disturbance  travels 
downstream.  The  growth  rate  for  the  n  =  4  mode  at  a  swirl  ratio  ofq  =  0.2  is  shown  in  figure  12. 
The  trends  are  similar  to  the  n  =  3  mode  presented  in  the  previous  figure  except  this  mode  attains 
its  peak  growth  rate  even  farther  downstream  at  z  =  160.  The  effect  of  increased  rotation  on  these 
higher  n  modes  is  very  similar  to  the  n  =  2  instability  and  therefore  are  omitted  here.  For  our  fixed 
Reynolds  number  of  Re  =  30000,  the  ra  =  1  mode  attains  the  largest  growth  rate  in  the  azimuthal 
wave  number  space.  The  absolute  maximum  growth  rate  is  expected  to  shift  to  larger  values  of  n  if 
Re  is  increased. 

For  streamwise  vortices,  in  addition  to  the  primary  unstable  mode  there  are  always  higher 
modes  of  instability  associated  with  each  individual  n .  This  is  irrespective  of  azimuthal  wave 
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number  considered.  As  the  Reynolds  number  is  increased,  the  number  of  these  higher  modes 
tends  to  increase  until  in  the  limit  of  Re  oo  there  are  an  infinite  number  of  unstable  modes. 
Fortunately,  the  growth  rates  of  these  higher  modes  are  much  lower  than  those  of  the  primary 
modes  and  therefore  are  less  important. 

The  growth  rate  behavior  of  the  first  higher  mode  for  the  n  =  1  instability  at  g  =  0.2  is  presented 
in  figure  13.  This  particular  mode  displays  a  double  hump  peak  growth  rate.  The  second  hump 
associated  with  the  higher  frequencies  tends  to  stabilize  much  faster  as  the  disturbance  moves 
downstream.  As  mentioned  above,  overall  the  growth  rate  is  significantly  smaller  than  the 
primary  mode  (figure  4).  The  corresponding  higher  mode  for  the  n  =  2  instability  is  shown  in 
figure  14.  For  the  remainder  of  this  report,  we  will  omit  the  discussion  on  the  higher  modes. 

5.1.1  Viscous  effects 

As  discussed  in  the  previous  sections,  the  Reynolds  number  does  not  play  a  dominant  role 
because  an  inviscid  (centrifugal)  mechanism  is  responsible  for  the  instabilities.  The  effect  of 
Reynolds  number  is  even  less  important  at  the  moderately  high  Re  of  most  experiments.  It  is 
instructive,  however,  to  show  a  sample  computation  where  viscous  effects  are  presented. 

The  inviscid  nature  of  the  helical  unstable  modes  is  shown  in  figure  15  where  we  have  plotted 
the  effect  of  Reynolds  number  on  the  n  =  1  disturbances  at  two  streamwise  locations.  For  Re  <  10^ 
the  instability  is  fully  damped.  As  the  Reynolds  number  is  increased  gradually,  the  primary 
mode  is  first  destabilized  and  then  after  a  period  of  rapid  growth,  its  amplification  factor  (Oi 
as3m[iptotes  eventually  to  a  constant  level  (inviscid  limit).  The  second  mode  becomes  destabilized 
once  the  Reynolds  number  is  increased  beyond  Re  =  4,000.  The  disturbance  frequency,  for  both 
modes  remains  fairly  independent  of  the  Reynolds  number  (not  shown).  The  behavior  of  the  n  =  2, 
and  3  modes  are  essentially  similar  to  that  of  the  n  =  1  perturbation,  except  that  the  high-Reynolds 
number  behavior  is  shifted  to  higher  Re5uiolds  number.  Based  on  these  observed  trends,  we  fixed 
the  Re5niolds  number  at  the  experimental  value  of  30,000  for  the  entire  computational  study. 

5.1.2  Compressibility  effects 

The  compressibility  effects  on  the  n  =  1  mode  for  q  =  0.2  are  presented  in  figures  16-18.  Three 
different  Mach  numbers  M  OO  —  !•  6,  2,  and  3  are  considered.  Results  for  M  00—2.  49  already  were 
shown  in  figure  3.  Increasing  the  Mach  number  has  a  stabilizing  influence.  At  higher  values  of 
Moo  there  is  a  reduction  in  the  growth  rate  of  the  disturbances.  For  example,  at  the  location  2  =  80  the 
maximum  growth  rate  for  Moo  =  3.0  shows  nearly  40%  drop  relative  to  the  peak  growth  rate  at  Moo  = 
1.6.  Overall,  this  Mach  number  effect  is  consistent  with  the  trends  observed  for  free  shear-layers 
without  swirl. 
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5.2  Linear  PSE 

Using  the  parabolized  stability  equations,  we  now  attempt  to  determine  spatial  evolution  of  the 
disturbance  field.  The  linear  PSE  provides  a  quick  way  of  mapping  the  underlying  vortex 
stability  characteristics  and  to  isolate  the  effects  of  mean  flow  changes  in  the  streamwise  direction 
(non-parallel  effects)  on  the  disturbances.  Once  again  the  reference  point  for  this  vortex  is  the 
axis3nnmetric  wake  where  q  =  0. 

The  variation  of  the  spatial  growth  rate  vs  frequency,  for  the  /i  =  1  mode  of  the  axisym- 
metric  wake,  at  several  streamwise  locations  is  presented  in  figure  19.  As  in  the  case  of  the  quasi- 
cylindrical  analysis  of  the  previous  section,  the  instability  maintains  a  nearly  constant 
amplification  rate  along  the  streamwise  direction.  Similar  curves  for  the  n  =  1  instability  of  the 
vortex  at  q  =  0.1,  0.2,  and  0.3  are  plotted  in  figures  20a-20c.  Initially,  increasing  the  swirl 
parameter  q  is  destabilizing  and  will  provide  higher  growth  rates  up  to  q  =  0.2.  Beyond  this  point, 
rotation  becomes  a  stabilizing  force  causing  a  lower  growth  rate  for  the  instability.  In  addition, 
the  growth  rate  diminishes  much  faster  along  the  streamwise  direction  at  higher  q's  (figure  20c). 
The  corresponding  plots  for  the  azimuthal  wave  numbers  n  =  2,  and  3  are  shown  in  figures  21a-21c 
and  22a“22c  respectively.  The  trends  observed  in  these  figures  are  identical  to  those  described  for 
the  n  =  1  instability.  Overall,  in  the  swirl  parameter  space,  the  absolute  maximum  growth  rate 
occurs  at  q  =  0.2  which  corresponds  to  the  values  obtained  from  an  analysis  of  Smart’s  et  al  (1994) 
experiment.  Both  n  =  l  and  n  =  2  instability  modes  have  comparable  growth  rates  at  this  q  with  n  = 
1  being  slightly  higher.  The  value  of  q  where  this  peak  amplification  occurs  depends  on  the  axial 
velocity  deficit  5,  which  was  chosen  to  be  0.35  to  match  the  experiment.  Obviously,  changing  S 
would  cause  q  to  shift  too.  Nevertheless,  the  trends  observed  so  far  will  re-emerge  once  a 
comparable  ratio  of  q/S  is  attained. 

The  eigenfunction  variations  across  the  vortex  core  for  the  n  =  1,  2,  and  3  modes  are  plotted  in 
figures  23a-23c.  The  eigenfunctions  were  determined  for  frequencies  corresponding  to  the 
respective  instability  maximum  growth  rates  -ajmax*  Furthermore,  they  are  normalized  with 
respect  to  the  maximum  axial  velocity  component,  These  helical  disturbances  are  confined 

to  the  vortex  core  and  lie  very  close  to  the  axis  of  the  vortex.  In  figure  23a,  the  n  =  1  mode  is  the  only 
mode  with  nonzero  fluctuations  on  the  vortex  axis.  This  particular  helical  mode,  usually  called  a 
flapping  mode,  plays  a  dominant  role  in  the  development  of  axis5nnmetric  free-shear  flows.  For 
the  0  perturbations,  the  temperature  fluctuation  dominates.  Figures  23a-23c  show  that  as  n  is 
increased,  the  radial  position  where  the  eigenfunctions  are  centered  also  increases.  In  this 
respect,  disturbances  with  relatively  high  azimuthal  wave  numbers  may  be  regarded  as  becoming 
“ring”  modes;  that  is,  the  perturbation’s  energy  is  concentrated  in  a  narrow  ring  of  fluid  (thus  the 
name  “ring”  mode)  situated  at  the  critical  layer  radial  distance  from  the  vortex  axis.  The 
pressure  fluctuation  is  not  plotted  in  figures  23a-23b  for  two  reasons.  First,  the  pressure 
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amplitude,  based  on  the  scale  shovm,  is  nearly  zero.  Second,  the  pressure  decays  very  smoothly 
and  has  a  nonoscillatory  behavior  as  r  becomes  large,  thus  indicating  that  there  is  no  energy 
radiation  to  farfield. 

5.2.1  Compressibility  effects 

The  variation  of  the  growth  rate  -a^  vs  frequency  (o  for  the  axisymmetric  wake  is  presented 
in  figure  24.  Results  were  obtained  for  five  distinct  supersonic  Mach  numbers.  The  fully- 
developed-wake  instability  tends  to  be  insensitive  to  the  freestream  Mach  number.  In  contrast,  the 
vortex  disturbance  field  is  quite  sensitive  to  Mach  number  variations.  Results  forn  =  1,  2,  and  3  at 
q  =  0.2,  2:  =  160  are  displayed  in  figures  25-27.  Compared  to  the  wake,  vortex  instabilities  go 
through  a  significant  reduction  in  their  growth  rate  when  the  freestream  Mach  number  is 
increased.  There  is  also  a  noticeable  reduction  in  the  range  of  frequency  where  the  instabilities 
are  present.  The  severity  of  the  compressibility  effects  is  slightly  more  pronounced  for  the  n  =  1 
mode  as  compared  to  the  other  two  azimuthal  wave  numbers.  This  point  is  made  more  clear  in 
figure  28  where  the  maximum  amplification  factor  max  plotted  vs  Mach  number.  Notice  at  M 
=  1.5  the  three  disturbances  have  comparable  growth  rates.  As  the  fi^eestream  M  is  increased,  the 
maximum  growth  rate  associated  with  the  n  -  1  disturbance  (the  curve  with  circles)  diminishes  at 
a  faster  rate  than  the  corresponding  curves  for  n  =  2,  and  3.  There  is  a  crossing  between  m  =  1  and  n 
=  2  that  takes  place  at  M  =  1.8  while  the  crossing  between  n  =  2  and  3  occurs  at  M  =  3.5.  As  indicated 
earlier,  the  wake  mode  is  virtually  unaffected  and  remains  nearly  constant  for  the  range  of  Mach 
number  considered.  Surprisingly,  figure  28  suggests  that  for  freestream  Mach  numbers  greater 
than  4  there  is  a  possibility  that  the  growth  rates  associated  vnth  vortex  instabilities  will  dip  below 
that  of  an  axis5mametric  wake. 

5.2.2  Non-parallel  effects 

The  linear  PSE  formulation  allows  one  to  determine  the  effects  of  streamwise  mean  flow 
variation  on  the  growth  of  the  disturbance  field.  Unlike  quasi-cylindrical  flow  (section  5.1), 
however,  there  is  no  unique  streamwise  wave  number  a  in  the  nonparallel  flow  analysis.  That 
is,  depending  on  the  path  taken  in  r-z  plane,  one  obtains  a  different  rate  of  amplification.  To 
circumvent  this  problem,  in  the  nonparallel  studies,  the  wave  number  can  be  designated  as  the 
“wave  number  based  on  a  particular  fluctuating  quantity  such  as  velocity,  temperature  or  mass- 
flow  rate.”  Usually  the  wave  number  is  computed  at  the  location  in  the  shear  layer  where  the  flow 
fluctuations  become  maximum.  For  example,  if  we  choose  to  compute  the  wave  number  a  at  the 
location  where  mass  fluctuation  peaks,  a  is  obtained  from  the  relation 
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where  r  is  the  radical  location  of  Accordingly,  the  growth  rate  then  is  defined  as 
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The  effect  of  the  nonparallel  mean  flow  on  the  /i  =  1  perturbation  with  q  =  0.2  and  ft>  =  1  is  shown  in 
figure  29.  We  have  chosen  four  different  fluctuating  quantities  namely  mass-flow  rate,  stream- 
wise  velocity,  temperature,  and  kinetic  energy  to  deduce  the  growth  rates.  In  figure  29  the  solid 
curve  represents  the  quasi-cylindrical  (parallel)  flow  analysis.  Independent  of  the  quantity 
considered,  the  computed  nonparallel  curves  all  follow  the  parallel  curve  and  are  virtually 
indistinguishable  from  it.  The  corresponding  plots  for  /i  =  2,  and  3  modes  are  shown  in  figures 
30-31.  For  these  azimuthal  wave  numbers  the  nonparallel  curves  show  slightly  smaller  growth 
rates  than  the  parallel  flow  curve.  Overall  the  differences  are  not  significant  which  suggests  that 
nonparallel  effects  do  not  alter  the  stability  characteristics  of  supersonic  vortices. 

5.3  Nonlinear  PSE 

The  linear  studies  provided  ample  information  to  choose  the  most  appropriate  test  cases.  The 
mean  flow  conditions  chosen  for  our  study  are  =  2.49,  Re  =  30,000,  Pr  =  0.76,  q  =  0.2  and  5  = 
0.35.  A  possible  scenario  for  laminar-turbulent  transition  is  through  interaction  of  the  vortex  with 
the  first  few  helical  modes.  The  initial  growth  of  these  disturbances  will  result  in  the  emergence 
of  organized  large-scale  structures.  Through  nonlinear  interactions  the  helical  modes  generate 
other  harmonics  which  lead  to  smaller  scale  structures  and  hence  transition  to  turbulence. 

For  the  present  study,  three  cases  were  computed.  For  the  first  case,  the  vortex  was  perturbed  by 
the  helical  mode  n  =  1  at  a  frequency  (o  =  1.0.  Because  the  Fourier  expansion  has  the  form 
ex^[in9-imojt),  we  denote  the  corresponding  mode  as  the  {m,n)  mode,  hence  the  initial 
disturbances  in  the  first  case  are  associated  with  the  (1,  1)  mode.  The  second  case  involves  the  =  2 
mode  which  was  introduced  at  a  frequency  of  1.3.  This  mode  is  designated  (1,  2)  with  co  =  1.3.  The 
last  case  considered  the  modal  interaction  between  n  =  l^n-2  and  the  vortex  at  co  =  1.3.  For  this 
case  the  initial  disturbances  are  the  modes  (1,  1)  and  (1,  2).  For  all  three  cases,  the  initial 
amplitude  for  the  axial  velocity  perturbation  was  taken  to  be  0.1%,  All  the  nonlinear  computations 
were  performed  with  expansions  to  M  =  5  and  N  =  5  with  M  and  N  defined  in  Eq.  4.24. 


27 


5.3.1  Case  a:  n=l  mode 

Figure  32  shows  the  streamwise  evolution  of  the  axial  velocity  fluctuation.  Throughout  the 
linear  stage,  the  (1,1)  mode  dominates.  The  growth  of  the  (1,1)  mode  follows  the  linear 
amplification  up  to  z  «  28;  further  downstream  the  growth  rate  departs  firom  the  linear  solution  and 
soon  nonlinear  saturation  takes  place.  In  the  nonhnear  stage  the  mean-flow  distortion  mode  (0,0) 
generated  by  nonlinear  interaction  overtakes  the  (1,1)  mode  and  rises  to  very  large  amplitudes 
(greater  than  7%)  before  it  too  saturates.  Prior  to  2  =  38,  the  higher  harmonics  do  not  attain  large 
enough  amphtudes  to  be  of  any  significance. 

The  various  components  of  the  mean-flow  distortion  mode  (radial,  tangential,  and  axial 
velocities  and  temperature)  are  plotted  in  figures  33a— 33d.  The  radial  velocity  (figure  33a)  is 
relatively  small  and  does  not  change  much  in  the  streamwise  direction.  The  tangential  velocity 
(figure  33b)  shows  moderate  growth.  On  the  other  hand  both  axial  velocity  and  temperature 
(figures  33c— 33d)  show  significant  levels  of  distortion  in  the  streamwise  direction,  in  particular 
during  the  saturation  stage.  The  axial  velocity  distortion  tends  to  reduce  the  wake  velocity  deficit. 
Meanwhile,  the  temperature  distortion  mode  reduces  the  temperature  excess  in  the  vortex  core. 
The  axial  velocity  contour  plots  in  the  r  -z  plane  are  shown  in  figure  34.  Towards  the  end  of  the 
linear  stage,  at  2  =  26,  the  instability  begins  to  modulate  the  vortex  core,  which  results  in  the 
appearance  of  a  large  scale  coherent  structure.  At  the  same  location,  the  reduction  in  the  wake 
velocity  deficit  starts  to  accelerate.  The  contour  plots  in  the  r  —  B  plane  at  2  =  1,  28,  and  38  are 
displayed  in  figures  35a-35c.  In  these  figures  the  solid  horizontal  line  is  the  branch  cut  in  B 
highlighting  B  =  0  line.  At  2  =  1,  the  vortex  is  fully  axisymmetric  and  the  contours  are  concentric. 
At  2  =  28,  the  evolution  of  the  disturbance  has  displaced  the  vortex  core  resulting  in  a  slight 
eccentricity  of  the  velocity  contours.  Near  the  end  of  saturation  stage  at  2  =  38,  the  core  is  displaced 
significantly  and  is  no  longer  axisymmetric. 

The  disappearance  of  the  axial  velocity  deficit  in  the  core  is  best  shown  in  figure  36  where  we 
have  plotted  the  iso-surface  of  velocity  at  a  constemt  value  of  0.74.  Recall  that  the  initial  starting 
velocity  has  a  35%  deficit.  Figure  36  shows  the  initial  twisting  and  shrinking  of  the  core  as  the 
disturbance  amplifies.  During  the  saturation  stage  the  wake  is  broken  into  little  islands  of  low 
axial  velocity  until  it  finally  disappears.  For  the  sake  of  completeness,  the  velocity  vectors  in  the 
cross  plane  at  z  =  1  and  39  are  presented  in  figures  37a  and  37b.  The  velocity  vectors  at  the  initial 
streamwise  station  show  the  vortex  in  its  axisymmetric  state  while  at  2  =  39  they  show  the  complex 
asymmetric  nature  of  the  velocity  field  before  breakdown  occurs. 

In  figure  38  the  contour  plot  for  the  temperature  field  in  r-z  plane  is  displayed.  Similar  to  the 
velocity,  the  temperature  field  also  clearly  shows  the  undulation  of  the  vortex  core  by  the 
disturbance  and  subsequent  generation  of  large-scale  coherent  structures.  Notice  the  vortex  starts 
with  a  temperature  excess  of  70%.  During  the  final  nonlinear  stage  the  core  temperature  starts  to 
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drop  significantly  in  the  streamwise  direction.  The  hotter  fluid  shows  up  in  smaller  regions  off  of 
the  centerline  in  the  middle  of  the  large  structures.  The  corresponding  contour  plots  in  the  r-0 
plane  at  2:  =  1,  28,  and  38  are  presented  in  figures  39a-39c.  Note  how  the  initial  axisymmetric 
temperature  field  becomes  eccentric  and  finally  develops  into  an  asymmetric  complex  distorted 
field.  At  the  2  =  38  location,  the  circular  core  has  been  deformed  into  a  strip  structure  which  is 
identified  by  the  contour  levels  F  and  E.  Once  again  the  breakdown  of  the  hot  vortex  core  is  best 
shown  by  the  iso-surface  of  temperature  displayed  in  figure  40  for  a  constant  level  of  1.54.  In  this 
figure  the  growth  of  the  instability  and  the  resulting  flapping  of  the  core  is  very  clear.  As  in  the 
case  of  the  velocity  field  (figure  36),  during  the  final  nonlinear  stage  the  core  breaks  down  into 
smaller  segments. 

Finally,  the  iso-surface  of  axial  vorticity  at  a  value  of  0.078  is  shown  in  figure  41.  Initially, 
the  axial  vorticity  resides  entirely  in  a  small  region  inside  of  the  core.  Due  to  the  l/r  behavior  of 
the  tangential  velocity,  the  flowfield  surrounding  the  core  is  basically  devoid  of  any  axial 
vorticity.  Figure  41  shows  that  axial  vorticity  follows  a  helical  path  in  accordance  with  the 
instability  wave  form.  Towards  the  beginning  of  the  saturation  stage,  the  shape  of  the  surface 
starts  to  deviate  from  a  clean,  well-defined  wave  form.  Due  to  nonlinear  interactions,  smaller 
scale  structures  begin  to  appear,  indicating  the  end  of  the  organized  core. 

5.3.2  Case  b:  n=2mode 

The  evolution  of  the  axial  velocity  fluctuation  for  this  helical  mode  is  shown  in  figure  42.  The 
forced  (1,2)  mode  grows  linearly  up  to  z  =  34  and  is  the  dominant  mode.  The  mean-flow  distortion 
mode  (0,0)  begins  with  a  much  smaller  amplitude  but  amplifies  at  a  higher  growth  rate.  Towards 
the  beginning  of  the  nonlinear  stage,  the  (0,0)  mode  overtakes  (1,2)  while  the  latter  mode  goes 
through  saturation.  The  mean  flow  distortion  mode  grows  to  a  very  large  amplitude  (more  than 
10%)  before  the  computation  was  stopped  due  to  lack  of  convergence.  The  higher  harmonics  never 
attain  a  large  enough  amplitude  to  be  of  any  significance. 

The  evolution  of  various  components  of  the  mean  flow  distortion  mode  are  presented  in 
figures  43a-43d.  As  in  the  case  of  the  n  =  1  perturbation,  the  axial  velocity  and  temperature  fields 
show  the  highest  levels  of  modification.  Unlike  n  =  1  mode,  however,  the  maximum  distortion 
here  occurs  off  the  vortex  centerline  at  a  radial  distance  of  r  =  0.25  and  0.6, 

The  axial  velocity  contour  plots  in  the  r-z  plane  are  displayed  in  figure  44.  Unlike  the  n  =  1 
disturbance,  the  vortex  core  remains  on  the  centerline  albeit  deformed  in  shape.  The  large  scale 
structures  inside  the  core  appear  at  2  =  38  and  grow  rapidly  downstream.  The  wake  velocity  deficit 
also  is  reduced  beyond  this  axial  location  but  not  as  severely  as  in  the  case  of  the  n  =  1  disturbance. 
The  velocity  contours  in  the  r-G  plane  at  2  =  30  and  50  are  shown  in  figures  45a  and  45b.  The 
contours  at  the  initial  station  z  =  1  are  identical  to  those  presented  in  figure  35a  and  therefore  are 
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not  repeated  here.  Figure  45a  shows  the  contours  towards  the  end  of  the  linear  growth.  As 
mentioned  above,  the  core  is  deformed  but  not  displaced  from  the  centerline  r  =  0.  Downstream, 
due  to  nonlinear  interaction,  the  core  is  severely  deformed  and  develops  into  an  s  shape.  In  other 
words,  the  vortex  core  is  being  deformed  into  a  two-lobed  structure.  This  pattern  becomes  typical 
for  disturbances  with  n  >  1.  That  is,  had  we  perturbed  the  vortex  with  n  =  3  mode,  the  core  would 
have  deformed  into  three  lobes  and  so  on.  Figure  46  shows  the  iso-surface  of  axial  velocity  at  a 
value  of  0.8.  The  appearance  of  the  disturbance,  its  subsequent  growth,  and  the  development  of  a 
two-lobed  core  structure  is  clearly  demonstrated  in  this  figure.  Notice  towards  the  end  of 
nonlinear  saturation,  the  surface  has  become  jagged  which  indicates  the  formation  of  smaller 
scale  structures. 

In  figures  47-49  the  corresponding  results  for  the  temperature  field  are  plotted.  The  trends 
just  described  for  the  axial  velocity  not  only  are  prevalent  for  the  temperature  field  but  actually 
show  up  in  a  more  pronounced  way.  Finally,  the  iso-surface  of  axial  vorticity  at  a  value  of  0.0795 
is  displayed  in  figure  50.  The  axial  vorticity  shows  the  presence  of  two  large  structures  that  grow 
in  the  downstream  direction.  The  structures  which  behave  as  two  distinct  intertwined  waves  start 
out  as  thin  vorticity  sheets.  As  the  disturbance  amplifies  downstream,  the  vorticity  sheets  grow  in 
size  and  become  thicker. 

5.3.3  Case  c:  Interaction  ofn  =  1,  and  2  Modes 

In  this  case  both  7i  =  1  and  2  modes  were  forced  at  a  m  =  1.3.  The  initial  forcing  level  was  kept 
to  0.1%  for  the  axial  velocity.  Figure  51  presents  the  downstream  evolution  of  various  Fourier 
modes.  The  two  forced  modes  (1,1)  and  (1,2)  amplify  in  the  linear  stage  with  nearly  identical 
growth  rates  and  remain  dominant.  As  the  nonlinear  stage  sets  in,  the  (1,1)  mode  saturates  at  a 
lower  amplitude  and  after  a  distance  of  six  core  radii  downstream,  its  amplitude  starts  to  drop  off. 
On  the  other  hand,  the  (1,2)  mode  saturates  and  levels  off  at  a  slightly  higher  amplitude.  Similar  to 
the  previous  cases,  the  (0,0)  mean  flow  distortion  mode  begins  at  a  fairly  small  amplitude. 
Towards  the  end  of  the  linear  stage  {z  =  28),  the  distortion  mode  overtakes  both  forced  modes  and 
reaches  very  high  amplitudes  (greater  than  10%).  Also  note  the  presence  and  importance  of  (0,1) 
and  (1,0)  modes  in  this  case.  Both  of  these  harmonics  attain  higher  amplitudes  than  the  (1,1)  forced 
mode  in  the  nonlinear  regime.  The  (1,0)  mode  represents  an  axisymmetric  traveling  wave. 

Figures  52a-52d  display  the  streamwise  evolution  of  the  mean  flow  distortion  mode. 
Relatively  speaking,  the  radial  and  tangential  velocities  remain  small.  As  before,  most  of  the 
distortion  is  related  to  the  axial  velocity  and  temperature  fields  with  temperature  modifications 
reaching  very  high  levels.  The  axial  velocity  contour  plots  in  the  r-z  plane  are  shown  in  figure  53. 
The  undulation  of  the  vortex  core  and  wake  deficit  reduction  begins  near  z  =  28.  Further 
downstream,  the  undulation  loses  its  cohesiveness  and  other  structures  start  to  form.  The  contour 
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plots  in  the  r-6  plane  at  two  streamwise  locations  0  =  30  and  40  are  presented  in  fibres  54a  and 
54b.  Note  that  by  2:  =  40,  the  vortex  core  has  evolved  into  a  complex  structure  that  neither  resembles 
the  71  =  1  mode  eccentric  shape  (figure  35c)  nor  does  it  have  the  s  shape  feature  of  the  n  =  2  mode 
(figure  45b).  The  above  point  is  made  more  clear  in  figure  55  where  we  have  plotted  the  iso-surface 
of  axial  velocity  at  a  value  of  0.74.  The  shape  of  the  surface  indicates  that  neither  of  the  two  forced 
modes  (1,1)  and  (1,2)  dominate  the  evolution  of  axial  velocity. 

The  temperature  contours  in  the  r-z  plane  is  shown  in  figure  56.  The  growth  of  instability 
waves  on  the  vortex  core  and  subsequent  formation  and  manifestation  of  large  scale  coherent 
structures  is  more  evident  in  this  figure  than  in  figure  53.  The  contour  plots  in  the  r-0  plane  are 
displayed  in  figures  57a  and  57b.  Atz  =  30,  the  hotter  fluid  in  the  core  has  been  displaced  off  the 
centerline  and  is  deformed  into  a  crescent  shape.  Atz-  40,  the  hotter  fluid  region  is  deformed  even 
more  and  resides  farther  away  from  the  centerline.  In  addition,  multiple  other  regions 
containing  warmer  fluid  have  appeared.  The  corresponding  iso-surface  for  the  temperature  is 
shown  in  figure  58.  Contrary  to  the  axial  velocity  iso-surface  (figure  55),  the  temperature  field 
shows  a  remarkable  resemblance  to  the  iso-surface  of  the  n  -  1  helical  mode  (figure  40).  The 
distortion  of  the  vortex  core  into  a  corkscrew  figure  is  apparent.  Also  apparent  is  the  breakdown  of 
the  core  in  the  nonlinear  stage  into  smaller  segments  of  hotter  fluid. 

In  figure  59  we  have  plotted  the  iso-surface  of  axial  vorticity.  In  the  linear  stage,  the  distortion 
of  the  vorticity  surface  initially  follows  a  path  similar  to  that  observed  for  the  n  =  1  disturbance 
(figure  41).  In  the  nonlinear  stage,  however,  the  interaction  between  the  two  forced  modes  has 
altered  the  surface  somewhat  with  the  smaller  scale  features  being  more  prominent. 

6.  Conclusions 

Quasi-cylindrical  linear  stability  analysis  was  used  to  thoroughly  map  the  stability 
characteristics  of  supersonic  streamwise  vortices.  Important  trends  such  as  viscous  and 
compressibility  effects  were  studied.  Based  on  linear  analysis,  the  unstable  disturbances  were 
found  to  be  inviscid  in  nature.  Except  for  relatively  low  Re3molds  numbers,  viscous  forces  were 
found  to  be  unimportant.  In  contrast,  compressibility  manifests  itself  as  a  stabilizing  force. 
Increasing  the  Mach  number  significantly  suppressed  the  growth  rate  of  the  helical  unstable 
modes.  The  reduction  was  more  pronounced  for  the  first  helical  mode  n  =  1. 

Non-parallel  and  nonlinear  effects  on  the  evolution  of  the  disturbance  field  were  studied 
using  linear  and  nonlinear  Parabolized  Stability  Equations  (PSE)  formulations.  Streamwise 
mean  flow  variation  was  determined  to  have  no  significant  effect  on  the  amplification  of  the 
instabilities.  For  the  nonlinear  analysis,  three  distinct  cases  namely  a)  71  =  1  mode,  b)  ti  =  2  mode, 
and  c)  71  =  1,  and  2  modal  interaction  were  considered.  In  all  three  cases,  nonlinear  growth  of  the 
disturbance  field  resulted  in  the  appearance  of  coherent  large  scale  motion  and  significant  mean 


31 


flow  distortion  in  the  axial  velocity  and  temperature  fields.  The  resulting  mean  flow 
modifications  tend  to  reduce  the  vortex  axial  velocity  deficit  and  temperature  excess  in  the  core. 
The  disappearance  of  axial  velocity  deficit  has  a  significant  implication  on  the  vortex  stability. 
Based  on  past  and  present  investigations,  the  presence  of  unstable  modes  in  the  vortex  core  is 
determined  almost  exclusively  by  the  ratio  of  the  swirl  parameter  q  to  the  velocity  deficit  (or 
excess)  5.  The  present  PSE  analysis  suggests  that  during  nonlinear  evolution  of  unstable 
disturbances,  S  is  reduced  while  q  remains  unaffected.  This  tends  to  increase  q  18  ratio.  The 
apparent  increase  in  the  swirl  level  provides  a  strong  stabilizing  force  which  will  revert  the  vortex 
core  back  towards  stability  and  a  laminar  state.  This  result  seems  to  support  conclusions  drawn 
by  Ragab  and  Sreedhar  (1995)  from  their  large-eddy  simulation  of  incompressible  streamwise 
vortices. 
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Table  I.  Flow  conditions  and  parameters  used  in  the  present  study. 


Free-stream  static  temperature, 

130°  K 

Free-stream  Reynolds  number,  Re  =  /  v,.. 

=  30,000 

Free-stream  Mach  number,  M„  =W^/a^ 

1.5  <  <  3.5 

Prandtl  number,  Pr  =  v„Cp/K„ 

0.76 

Swirl  parameter,  q 

0<q<0A 

Axial  velocity  deficit,  S 

0.35 

35 


q=0.2 ,  M„=2.49 


u 

Fig.  Ic.  Radial  velocity. 


Fig.  Id.  Total  temperature 


Fig.  3.  Variation  of  the  vortex  instability  growth  rate  vs.  frequency  with  =  2.49,  Pr  =  0.76, 
and  Re  =  3x10^. 


Fig.  4.  Variation  of  the  vortex  instabihty  growth  rate  vs.  frequency  for  Af„  =  2.49,  Pr  =  0.76,  and 

Re  = 3x101 


40 


)f  the  vortex  instability  growth  rate  vs.  frequency  for  =  2.49,  Pr  =  0.76,  and 


j-1 _ 


vs.  freauencv  for  =  2.- 


SjTT 


n=2 ,  q=0.1 


;he  vortex  instability  growth  rate 


®r 


Growth  rate  behavior  of  higher  unstable  mode  vs.  frequency  for  =  2.49,  Pr  =  0.76, 

andi?e  =  3x10"^. 


Growth  rate  behavior  of  higher  unstable  mode  vs.  frequency  for  =  2.49,  Pr  =  0.76, 
andiJe  =  3x10"^. 


5.  The  effect  of  increasing  the  Reynolds  number  on  the  growth  rate  of  the  n  =  1  unstable 
helical  mode. 


Variation  of  the  vortex  instability  growth  rate  vs.  frequency  with  Pr  =  0.76  and 
Re  =  3x104. 


Fig.  17.  Variation  of  the  vortex  instability  growth  rate  vs.  frequency  with  Pr  =  0.76  and 

Re  =  3x104. 


Fig.  18.  Variation  of  the  vortex  instability  growth  rate  vs.  frequency  with  Pr  =  0.76  and 

Re  =  3x104. 
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Variation  of  the  disturbance  spatial  growth  rate  vs.  frequency  for  the  axisymmetric  wake 


(a)  Low  swirl 

Variation  of  the  spatial  growth  rate  of  the  first  helical  mode  vs.  frequency  for  different 
values  of  swirl  parameter  q. 


Fig.  21. 


0) 

(a)  Low  swirl 

Variation  of  the  spatial  growth  rate  of  the  second  helical  mode  vs.  frequency  for  different 
values  of  swirl  parameter  q. 


Fig.  21b.  Moderate  swirl 
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Amplitude 


Radial  distribution  of  the  eigenfunctions  for  the  first  three  helical  unstable  modes  of  the 
vortex. 
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M=2.49,  Re=30000,  q=0.2,  z=160,  n=2,  co=1.3 
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The  effect  of  increasing  the  Mach  number  on  the  growth  rate  of  the  first  helical 
instability  mode. 


The  effect  of  increasing  the  Mach  number  on  the  growth  rate  of  the  third  helical 
instability  mode. 


(a)  Radial  velocity 

Distribution  of  various  components  of  meanflow  distortion  (0,0)  mode  in  the  radial 
direction. 


Fig.  33b.  Tangential  velocity. 
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M=2.49,  Re=30000,  q=0.2,  n=1,o)=1.0 


(a)  Initial  field. 

Fig.  37.  Cross  plane  velocity  vectors  showing  the  initial  axisymmetric  flowfield  and  the  end  of 
nonlinear  saturation  stage. 


Fig.  37b.  Nonlinear  saturation  stage. 
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Temperature  contour  plots  in  the  r-z  plane  showing  the  appearance  of  large  scale 
coherent  structures  in  the  nonlinear  stage. 


fso-Surface  of  axial  vorticity,  OMEGA_z  =  0.078 


Fig.  41.  Iso-surface  of  axial  vorticity  showing  the  evolution  of  vorticity  field  in  the  vortex  core. 
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Evolution  of  axial  velocity  fluctuations  with  axial  distance  showing  the  dominant  helical 
(1,2)  mode  along  with  other  Fourier  harmonic  waves. 


M=2.49,  Re=30000,  q=0.2,  n=2,  co=1.3 
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M=2.49,  Re=30000,  q=0.2,  n=2,  oo=1.3 
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Fig.  47. 
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Temperature  contour  plots  in  the  r-z  plane  showing  the  appearance  of  large  scale 
coherent  structures  in  the  nonlinear  stage. 
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(a)  z  =  30 

Fig.  48.  Temperature  contour  plots  in  the  r-0  plane  showing  the  loss  of  axisymmetry  as  the 
instability  evolves  downstream. 
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M=2.49,  Re=30000,  q=0.2,  n=1&2,  co=1.3 


Evolution  of  axial  velocity  fluctuations  with  axial  distance  showing  the  dominant  helical 
(1,1)  and  (1,2)  modes  along  with  other  Fourier  harmonic  waves. 


(a)  Radial  velocity 

Distribution  of  various  components  of  meanflow  distortion  (0,0)  mode  in  the  radial 
direction. 
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Fig.  53.  Axial  velocity  contour  plots  in  the  r-z  plane  showing  the  evolution  of  vortex  core  in  the 
streamwise  direction. 
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Fig.  54b.  z  =  40 


Temperature  contour  plots  in  the  r-z  plane  showing  the  appearance  of  large  scale 
coherent  structures  in  the  nonlinear  stage. 


Iso-Surface  of  axial  vorticity,  OMEGA_z=0.078 
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Fig.  59. 


Iso-surface  of  axial  vorticity  showing  the  evolution  of  vorticity  field  in  the  vortex  core. 


